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ABSTRACT 

In order for the NIF (National Ignition Facility) to achieve its most challenging 
goal of fusion, it is paramount that we explore new laser matter interaction regimes and 
control mechanisms. This thesis explores a novel] interaction regime namely, the 
enhancement of laser absorption when an overdense plasma is irradiated with ultra- 
intense laser light. Self-generated magnetic fields of up to 10° gauss have been predicted 
in overdense plasmas irradiated with ultra-intense laser light pulses. We find that laser 
light absorption can be significantly enhanced by the oscillation of electrons across these 
fields. There is then a very large Lorentz force, which can strongly accelerate electrons 
into the plasma and can lead to generation of harmonics in the reflected light. We also 
show that large magnetic fields can be significantly amplified to even larger values by the 
pressure of the light pulse. Potential applications of this interaction regime include the 
Fast Igniter approach to ICF (Inertial Confinement Fusion), the creation of anti-matter, 
and the development of laser radiography to name a few. Zohar, a 2-2 D PIC (Particle In 
Cell) code was used to run our simulations. Zohar assumes a collisionless plasma, solves 
Maxwell’s equations in 3D, and accounts for the relativistic effects which can occur in 
the high-energy regime that we consider. The high number of particles simulated 
required significant computing power, which was obtained by remote access to the LC 
(Livermore Computing) East computer, which provided 8 processors and 8 GB of shared 
memory at 440 MHz. “Exceed” software was required to graphically display our results 
on a PC. This thesis has been a collaborative effort between NPS (Naval Postgraduate 


School) and LLNL (Lawrence Livermore National Laboratories), to the benefit of DoD. 
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I. INTRODUCTION 


The study of the interaction of intense laser light with plasmas serves as an 
excellent introduction to the field of plasma physics [Ref. 1]. Both the linear and 
nonlinear theories of plasma waves, instabilities and wave-particle interactions are 
important for understanding laser plasma coupling. Numerous plasma effects have been 
observed in laser plasma experiments, and many challenging problems remain to be 


understood. 


A. BASIC PLASMA CONCEPTS 


Plasma is a medium containing many charged particles governed by 
electromagnetic forces. Plasma is a quasineutral gas of charged particles, which possess 
collective behavior. Plasma designates the fourth state of matter, which exists in the sun, 
stars, and space, etc; moreover, most of the Universe is made up of plasma. The plasma 
medium is described macroscopically (on a large scale) by its temperature and density, 
and changes in the plasma are calculated by using conservation equations such as 
conservation of energy, momentum and mass. On a microscopic (small) scale, a 
statistical description using probabilities is required to describe a plasma in the 
calculation of the positions and velocities of all particles. Due to their mutual collisions, 
the plasma charged particles emit radiation (electromagnetic waves). Moreover, many 
different waves can be created in a plasma medium. The stabilities and instabilities of 
these waves play a crucial role in plasma systems. Plasma plays an important role in 
explaining our universe and, as we shall see, scientists believe that the understanding of 


plasma is the key issue for solving the energy problems of mankind. 


A plasma is basically a system of N charges that are coupled to one another via 
their self-consistent electric and magnetic fields. Following the evolution of these N 
charges in principal would require us to solve 6*N coupled equations which Roel be 
tedious at best. Fortunately, great simplification is possible if we focus our attention on 
the collisionless plasma behavior. 
Let us assume a plasma space in which we decompose the electric field into two 
fields (E, and E>) that have distinct spatial scales. The field E; has spatial variations on a 


scale length much less than the electron Debye length 2,,. The field of an individual 


charge is shielded out by the response of the surrounding charges over this length. E; 
represents the rapidly fluctuating microfield due to multiple and random encounters 
(collisions) among the discrete charges. In contrast, the field E2 varies on a space scale 


greater than or comparable to 4,,. E2 represents the forces due to the average or 


“collective” motion of the charges. 

We thus have a natural separation into collisional and collective behavior. The 
collisional behavior becomes negligible when the number of electrons in the sphere with 
a radius equal to the electron Debye length becomes very large. 


1. Concept of Temperature 


In order to better understand the complex nature of plasma behavior, it is 
necessary to expand our most basic notion of temperature [Ref. 2]. We note that a gas in 
thermal equilibrium has particles of many different velocities and the most probable 
distribution of these velocities is known as the Maxwellian distribution. The average 
kinetic energy (Eaye) in three dimensions is: 


Bees 3/2 KT (1) 


where K is Boltzmann’s constant, K = 1.38 x 10° J/°K and T is the temperature in 
degrees Kelvin. The general result here is that Eaye is “% KT per degree of freedom. 

It is customary in plasma physics to give temperatures in units of energy. To avoid 
confusion on the number of dimensions involved, it is not Eye but the energy 
corresponding to KT that is used to denote the temperature. For KT = 1 eV = 1.6x1 or J, 


we have T = 1.6x107'9/1.38x10™ = 11,600, thus the conversion factor is 


leV = 11,600° K. 

It is interesting that plasma can have several temperatures at the same time. It 
often happens that the ions and electrons have separate Maxwellian distributions with 
different T; and T, (ion and electron temperatures respectively). This comes about 
because the energy exchange between colliding electrons is larger than that between an 
ion and electron. Since colliding electrons are the same size and mass, appreciable 
kinetic energy exchanges can result from an electron-electron collision. Now, when an 
electron-ion collision occurs, we recall that the electron is only 1/3700" the mass of an 
ion and therefore a finite amount of time (on the order of hundreds of microseconds let’s 
say) is required for a couple thousand electrons to collide with an individual ion to result 
in an appreciable kinetic energy exchange. Bearing this in mind, we can see how each 
species can be in its own thermal equilibrium but the plasma may not last long enough for 
the two temperatures to equalize (since many of the interesting effects of the experiments 
are observed in nanoseconds). 

When there is a magnetic field B, even a single species, say ions, can have two 
temperatures. This is because the forces acting on ions along B are different from those 


acting perpendicular to B (due to the Lorentz force). The components of velocity 


perpendicular to B and parallel to B may then belong to different Maxwellian 
distributions with temperatures Tparatte: ANd Tperpendicular- 

Interestingly enough, a high temperature does not necessarily mean a lot of heat 
[Ref. 2]. For example, the electron temperature inside a fluorescent light bulb is about 
20,000° K, but it is cool to the touch. This is because heat capacity has to be considered. 
The heat capacity of the electrons in the fluorescent light bulb is minimal because the 
electron density is much less than that of a gas at atmospheric pressure and the amount of 
heat transferred to the wall by electrons striking it at their thermal velocities is minimal. 


2. Plasma Instabilities 


Plasmas behave in a collective way. Different types of waves are created in the 
plasma due to the collective behavior of the electrons and the ions [Ref. 2]. The types of 
waves that can be created in plasma are numerous and we shall mention just a few of the 
most fundamental ones to the study of plasma behavior. 

Conditions of collective motion can be stable or unstable. For example, suppose 
a small displacement takes place in plasma confined by a magnetic field. If the 
collective motion behaves in a stable way, the displacement will be restored to the 
original motion. However, if the collective motion is unstable, any small displacement 
increases rapidly resulting in the collective motion phenomenon breaking up. 

Numerous instabilities can be excited by an intense electromagnetic wave in a 

plasma. Many possibilities exist even without dc magnetic fields [Ref. 1]. These 
various instabilities can be thought of as the resonant coupling into two other waves. 
Various possibilities and their names are shown below. 


= Ion acoustic decay instability a, > @ poo gota = Ne 


« Raman instability OPO, +O, rao Ys Nor 

« Brillouin instability DO, > O,.+@, ans No 

= 2@,, instability DO) > OO, +@,, atn= “Ne 
Here, @, is the incident light frequency, @,, is the frequency of a scattered light wave 
@,. is the frequency of a high frequency electron plasma wave, @, is the so-called ion 


sound frequency and n,, is the critical density in the plasma at which an incident light 
wave is reflected. 


3. Self Generated Magnetic Fields in Plasmas 


Experiments and calculations have shown the existence of self-generated 
magnetic fields [Ref. 3]. These magnetic fields can be generated in many different ways 
[Ref. 4]. Let’s consider the simplest way. Consider a target irradiated by high intensity 
laser light. There will be both a density and temperature gradient in the resulting plasma, 
as indicated in Figure 1. When the cross product between the density gradient (Vn) and 


the temperature gradient (VT) is nonzero, a magnetic field is generated. Faraday’s law 


shows that B is nonzero (i.e., a magnetic field grows up) whenever there is a curl of the 
electric field. In an expanding plasma, the electric field is proportional to the gradient of 
the pressure (VnKT) divided by the density. One then finds that a magnetic field is 


generated whenever Vn x VT #0. 


Target 


VT 
Incident | 
base a rr 
eit) 
—> 
Vn 


Figure 1. A schematic illustrating crossed density and temperature gradients 
generated when a target is irradiated. 


B. COMPUTER SIMULATION OF PLASMA USING PARTICLE CODES 


A theoretical description is often used where the plasma is treated as two charged 
fluids. We will use a complementary numerical description of plasma behavior using 
particle codes [Ref. 1,5,6]. Computer simulation of plasma using particle codes is a very 
direct and powerful approach, particularly for investigating kinetic and/or nonlinear 
effects. Since we are investigating strongly nonlinear plasma effects, we utilize the 
power of the computer to simulate the behavior of the charged particles of plasma. 

The approach is as follows: numerically follow the motion of a large collection of 
charges in their self-consistent electric and magnetic fields. The basic cycle is illustrated 
in Figure (2). 

From the positions and velocities at any given time, we compute the charge and 
current densities on a spatial grid that is fine enough to resolve the collective behavior of 
the charges. Using the computed charge and current densities, then compute the self- 


consistent electric and magnetic fields via Maxwell’s equations. Next, we use these 


fields in the equations of motion to advance the positions and velocities of the charges. 
Finally, continue around this basic cycle with a time step sufficiently small to resolve the 


highest frequency in the problem (which is often the electron plasma frequency @,,, ). 


Newton’s Grid 
laws assignment 


Maxwell’s 
equations 


Figure 2. The Basic Cycle of a Particle Simulation Code, From Ref.[1]. 

While plasma actually consists of a huge number of particles, particle simulations 
are limited by the capabilities of computers to 10* to 10° particles, and a single simulation 
particle must therefore vepnieseti a large number of actual particles. Since we seek to 
describe the collective behavior of the plasma, we can follow behavior on characteristic 
length £, which is larger than the Debye length /,, , to simplify our analysis and not 
consider the fine-scale fluctuating microfields associated with particle discreteness. This 
simplification allows us to represent the fields on a discrete grid in space, which is fine 
enough to resolve the collective behavior and sufficiently coarse to ignore the 


microfields. A mesh size 6 ~ /,, is usually chosen. 


Alternatively, this spatial averaging process can be viewed as using “finite-size” 


particles of size 5 ~ A,,. The behavior of a group of such particles is the same (with 


minor modifications) as the behavior of point charges for scale lengths ¢>> 6, but 
fluctuations with scale lengths £<< 6 are suppressed. Hence, we achieve with a “trick” 
what nature does with the use of an enormous number of particles; 1.e., smooth over the 


fine scale length microfields. 


1. Basic Ingredients of a Particle Code 

All particle codes share some basic ingredients, which we shall illustrate. We 
treat the ae as a fixed, neutralizing background, assume no imposed magnetic fields, 
and consider only electrostatic fields. In this electrostatic limit, the magnetic field 
generated by plasma currents is negligible, and Maxwell’s equations reduce to Poisson’s 
equation, which for simplicity we look at in only one dimension: 

E)V'$ = £09" /dx" =-e(n, —N,) (2) 

In this equation, €9 is the permittivity constant of free space, o is the electric field 


potential, e is the electron charge, n, and n, being the electron and ion charge densities 


respectively. Variations are allowed in only one direction, and periodic boundary 
conditions are adopted. 

Our first task is to compute the charge density. Our system extends from 0 to L, 
as Shown in Figure (3). We divide this system into NC cells and for convenience take the 
cell size (5) to be unity (i1.e., L= NC). The grid points are identified as the integer values 
of the position, augmented by one so that the counting begins with one. Note that the 
(NC + 1)" grid point is then identical with the first, due to the assumed periodic boundary 
conditions. Given the position of the charge, many schemes can be used to assign it to 


the spatial grid. For example, we could just assign the charge to its nearest grid point 


location. A better scheme is to share the charge between its two nearest grid points. (see 


Figure 3) 





cell 1 celli +] 


Figure 3. The charge-sharing scheme for finite sized particles, From Ref.[1]. 


For a charge located a distance A x to the right of the i grid point, we then have: 
Ap(i) = q(1 - Ax) (3) 
Ap(i +1) =q Ax , 
where Ap is the increment to the charge density. 
Having assigned the charges and determined the charge density on the spatial 
grid, we next determine the self-consistent electric field using Poisson’s equation 
OE/dx = pleo (4) 
The simplest approach is to finite-difference Poisson’s equation, 


E(i +1) = E@) + 6|(oi+1)+ o@)/2Ki/e.) (5) 


An alternative approach is to Fourier transform the charge density, use Poisson’s 
equation to find E, (ikE, = p,/€), and invert the transform to determine the electric field 
on the spatial grid. 
The last step in the basic cycle is to use the electric field to move the particles. 
The electrical force is assigned from the grid to the individual particles using the same 
scheme chosen to assign the charges to the grid. For example, considering a particle a 
distance Ax to the right of the i" grid point and using linear interpolation, we obtain for 
the force F on the particle 
F = qE(i){l—Ax]+ gE(i+1)Ax (6) 
In this code, Poisson’s equations are only solved for the electric field and the magnetic 
field is assumed zero. Therefore, the velocity (v) and position (x) of each particle can be 
advanced At in time using a “leap frog’’ algorithm; i-e., 
Veena yee ye” AF (7) 
x ay? py Pll? Ay (8) 
The superscripts denote time step. By defining x and v one-half time step apart, we 
achieve second-order accuracy in the time step. In the initial conditions, x and v are 
defined at the same time (t=0), but it is straightforward to then displace the velocities 
backward in time using the force at t=0. 
The plasma evolution is computed by simply continuing around the basic cycle 
using a time step small enough to resolve the characteristic oscillations of the plasma. In 


this electrostatic limit, the highest frequency oscillation has a frequency near @pe, the 


electron plasma frequency. Hence, a time step of about 0 20%, is commonly used 


Of course, it is also necessary to resolve the scale lengths characteristic of the collective 
behavior. Hence a grid size of about Age is commonly used. It is interesting to note that a 
numerical instability would be introduced if the grid size were chosen to be too many Age. 
This instability is due to aliasing, which arises from the fact that on the grid a disturbance 
with wave number k cannot be distinguished from spurious ones with wave numbers k + 


22n/6 where 6 is the grid spacing and n is an integer [Ref. 1,19]. 


2. PIC Code Programming Considerations 


Random particle initialization is the simplest and most commonly used method. 
However, it introduces noise in the computations which can mask interesting physical 
phenomena, particularly when we wish to study the low amplitude linear regime. In such 
cases the programs utilize what is referred to as a “quiet start” which avoids this 
background noise introduced by the randomness of the numbers. The program that 
simulates the plasma particle behavior can remain in this “quiet” state for only a finite 
time, which must be long enough to allow observation of the low-amplitude phenomena 
being studied. In such cases, staggering the particles, so that they do not form exact 
beams, is an effective way to delay the onset of the noisy behavior. 

Even with this simplest of particle code, many instructive and interesting plasma 
problems can be examined. As one example, let’s use ES1, a flexible and optimized 1d 
electrostatic PIC code to look at the strongly non-linear behavior of a two-stream 


instability [Ref. 5]. 


In this simple case, half the electrons are given an initial drift velocity in one 
direction and half in the other direction. This counter streaming motion is a source of 
free energy, which drives unstable a variation in density and potential. 

Figure 4a shows the electron phase space (velocity versus position for all the 
electrons), electric field and the potential at an early time. Note the counter streaming 
plasmas. In Figures 4b and Figure 4c, the effect of the instability on the electric field and 
the potential in the plasma can be seen. Finally, the growing electric field becomes so 
large that vortices form in phase space, as the electrons are alternatively slowed down 
and speeded by the large electric field, as shown in Figure 4c. Figure 4d shows time 
histories of the kinetic energy, the electric field energy and the total energy of the system. 
Note the exponential growth of the electric field energy at early times and the subsequent 
saturation. As the electric field energy increases, the kinetic energy decreases. The total 
energy varies very little compared to these changes, showing that the code is conserving 


energy with sufficient accuracy. 
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Figure 4c Figure 4d 





Figures 4a-c show V, Phase Space , Electric Field, and Potential respectively 
versus x. Figure 4 d shows’ Kinetic, Electric Field, and Total energies respectively versus 
time. 

Zohar represents a significant improvement to ES1 in PIC modeling of plasma 
behavior. This code solves all of Maxwell’s equations in order to include light waves and 
uses equations of motion that are relativistically correct. Zohar (which means brilliance 
or radiance in Hebrew) is the 2 % D PIC code used at Lawrence Livermore Labs to model 
the interaction of intense laser light with plasmas and for some electron beam-plasma 
studies[Ref. 7]. 2 % D means that variations are followed in two directions, and that al] 
three velocities are included in the code calculations. We note that state of the art, three- 
dimensional PIC codes, which run on parallel machines, now exist [Ref. 8,9]. The 


simulations that follow are results of Zohar runs. 


II. ULTRA-INTENSE REGIME OF LASER MATTER INTERACTIONS 


The recent development of Ultra Short Ultra-Intense (US-U]) laser pulses has 
allowed for the generation of laser intensities that are orders of magnitude over what was 
possible with previous laser technology [Ref.10]. The ultrashort regime includes laser 
light with pulse lengths t, of about 10 femtoseconds (fs) up to 10 picoseconds (ps). The 
ultra-intense regime covers intensities starting at about 10'” W/cm’. The ultra-intense 
regime is usually defined as the point at which the normalized amplitude of the electric 
field of the laser light is greater than or equal to 1, i.e., eE,/m.a@ oc = 1 [Ref.11]. Here, 
Eis the electric field amplitude of the laser light, m, is electron rest mass, @p 1s the 
frequency of the laser light, and c is light speed. Even as late as ten years ago, most 
lasers could not produce light that could approach this electric field strength, with the 
exception of a few CO, lasers built in the 1970’s. Now, there are a large number of ultra- 
intense lasers all over the world. What made this possible was the application of chirped 
pulse amplification (CPA) [Ref. 10], originally developed in microwave technology, to 
laser systems. This breakthrough in laser hardware has spawned considerable interest in 
US-UI laser-matter interactions. Never before has it been possible to deposit so much 
laser energy in such a short amount of time, in such a tiny volume on the front of a solid 
target. 

Chirped Pulse Amplification is critical because laser pulses of extremely high 
power density (gigawatts/centimeter’) can severely damage optical components such as 
amplifiers, lenses, and mirrors. CPA technology allows very short laser pulses to 


generate extremely high peak powers by stretching the initially low-energy laser pulse 


more than 10,000 times its duration before amplifying it and recompressing the pulse 
back to near the original duration after amplification. Damage to the laser optics is 


avoided by passing the long low power pulse through the amplifier (see Figure 5). 
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Figure 5. The concept of chirped-pulse amplification and other new technologies 


enabling the production of ultra-intense pulses, From Ref.[10]. 


Previously, the intense CO2 lasers were nanoseconds long, which creates large 
preformed plasmas that the laser interacted with. There was then a considerable amount 
of absorption in the underdense part of the plasma, as well as enough time for the laser to 
filament on its way to the critical surface. This, in turn, led to a large fraction (> 10%) of 
the laser energy converted to hot electrons, which were thought to be created by various 
Jaser-plasma interactions. However, US-UI laser pulses can have pulse lengths of the 


order of picoseconds and shorter. This means that they have the ability to both create a 


dense plasma (sometimes several times the critical density and up to solid densities) and 
interact with it, if proper care is taken to eliminate the pre-pulse. 

The interaction of an intense laser with very dense plasma is interesting for 
several reasons. First, since the laser light is absorbed most efficiently near the critical 
surface, the material can potentially absorb a large fraction of the laser light. Laser light 
transfers energy to the plasma via the electrons. Electrons readily oscillate in the electric 


field of the laser. This motion is referred to as the quiver velocity, and is defined as 


Pose |MpC =YVq,,/€ = CE, /m,c@, = 12. /1.3x 10% (9) 


Here, Posc (Vosc) 18 the transverse quiver momentum (velocity) of an electron in the laser 
field with an electric field amplitude of Eo, I the laser intensity (in W/cm”), Ay is the laser 
wavelength in microns, and y is the relativistic factor defined below. The energy 
associated with this momentum can easily reach MeV levels for electrons in the presence 


2 Therefore, relativistic effects 


of lasers with (normalized) intensities 10’? W-um?/em 
become important for these interactions, and are included in calculations concerning US- 
UI laser-matter interactions. Note that the quiver velocity scales in terms of IA”. This is 


a common scaling, and so it is useful to introduce this as the normalized intensity. In this 


way, many results discussed can easily be applied to any wavelength laser desired. 
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A. SIMPLE MODELS FOR REALISTIC LASERS AND SOLIDS 


An ideal model for ultra-intense pulses interacting with solids would require that a 
perfect Gaussian laser pulse, with a temporal pulse length (FWHM) of to and a peak 
intensity of J, strike a solid with a step-like density interface [Ref.11]. When the laser 
intensity exceeds about 10'* W/cm”, ionization of the electrons in the target occurs, and a 
plasma is formed on the surface of the solid. Typical ion densities are nj ~ 107? cm®, 
and electron densities are ne = Z nj. For US-UI pulses, this happens quite early in the 
interaction, and so most of the pulse interacts with a plasma, as opposed to an un-ionized 


solid. The relativistically correct dispersion relation (for circular polarization) that 


governs the propagation of light waves in the plasma is given by 


ke =,Jo2,/y- 0% (10) 


where k(q@, ) is the wavenumber (frequency) of the laser light, 


y = V1+ p2./(mey (11) 


is the relativistic factor, and 


Op, = 42n,e’ |m, (12) 


is the plasma frequency. 

First, consider the nonrelativistic case where y ~ 1, which is valid up to 
(normalized) intensities of about 5 x 10'” W-m7/cm”. In this case, k approaches zero 
When @pe= Yo. The laser no longer propagates, but is reflected when this condition is 
satisfied. The density at which this occurs is called the critical density and is given by 


yell 0 ere (13) 


where A, is the laser wavelength in micrometers. The laser light does penetrate slightly 
past critical, but falls off exponentially in a distance determined by the plasma density. 
This distance is known as the collisionless skin depth and is given by c/@p- in the 
nonrelativistic limit. For example, for a solid target of carbon (Z = 6), the skin depth is 
about 6 nm (nanometer = 10° m). 

For the relativistic case (y >1), we see that the effective density at which the light 
will be reflected from the plasma has increased by the factor y. Thus, in the relativistic 
case laser light can propagate even further into the overdense plasma due to the mass of 
the electrons increasing as they become relativistic. This increase in inertia makes it 
more difficult to generate the oscillating current in the plasma needed to reflect the laser 
light. 

A second interesting effect evident in US-UI laser matter interactions is that the 
laser pressure can easily be larger than the plasma pressure, even for extremely dense 
plasmas. The pressure associated with a laser with intensity J (in W /cm’) is 

P, = 3301/10'* Mbar (14) 

Efficient coupling of laser light to electrons has been exploited in two ways to 
date. First, there has been considerable effort at trying to use these interactions to 
produce high brightness X-ray lasers [Ref. 11]. Second, the possibility of generating 
large numbers of energetic electrons has been proposed as a novel way to deliver a large 
amount of energy to a hot spot in a compressed core of deuterium and tritium in order to 
ignite a fusion pellet. This is commonly known as the Fast Igniter fusion scheme 

[Ref. 12]. Central to both of these applications is the question of efficiency. How 


efficient is the coupling of laser light to the production of hot electrons? How does it 
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scale with laser intensity, wavelength, f-number, or spot size? How does the efficiency 
scale with the plasma parameters such as temperature, density, and scale length? [Ref.13] 
Computer simulations using Particle-In-Cell codes were used to provide an early 


glimpse into the ultra-intense regime. As a simple example, consider a two-dimensional 
simulation of a focused light beam with intensity lA, =10°W - um? /cm? incident onto 


an overdense plasma. Over a third of the light was found to be absorbed into very 
energetic electrons as shown by Figure 6a [Ref. 13]. These electrons had an effective 
temperature of about 1 MeV at this intensity. As shown in Figure 6b, the light pressure 
punched a hole in the plasma, causing it to recede at a significant fraction of the velocity 
of light. Both of these effects have now been observed. Very large self-generated 


magnetic fields were found in these simulations (up to 10° Gauss!). 
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Figure 6a, From Ref.[13]. A plot of the energy distribution of escaping electrons. 
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Figure 6b, From Ref.[{13]. A plot of electrons in space showing a hole being punched 


into the plasma. 
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Il. EFFECT OF LARGE SELF-GENERATED MAGNETIC FIELDS ON 
ULTRA INTENSE LASER PLASMA INTERACTIONS 
The focus of this thesis is to show how large self-generated magnetic fields affect 
the interaction of ultra-intense laser pulses with plasmas as found in the early simulation 
studies just discussed, the magnitude of the self-generated magnetic (Ba) field can be 
estimated as B,,/3, where Bos is the amplitude of the oscillating magnetic field of the 
short, intense laser pulse. In other words, about 10% of the energy of the oscillating 
magnetic field was “rectified”. When targets are irradiated with 1.06um laser light with 
intensities ranging from 10’? W/cm’ to 107! W/cm’, Bac can become very large, with 
typical values of 30-900 MG as shown in Table I. To put such large magnetic fields into 
perspective, we note the following: 
- The approximate magnetic field at the earths surface 
~ 10° T = 1 Gauss. 
- The approximate magnetic field near a big electromagnet 
~ 1.5 T= 1.5 x 10* Gauss. 
- The approximate magnetic field near the surface of a neutron star 


~ 10° T = 10? Gauss [Ref.14]. 
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Die ae, 





Table 1. Typical estimated values of Bg, and cyclotron frequency in the ultra- 
intense laser regime. Here @ce= eBa-/mcy and y = (1 + 1X, /2.8 x10'8 Pe 


Some other characteristic parameters are also shown in Table I. The electron 
motion in such intense laser light can become relativistic. The relativistic factor y is 


estimated as: 


y = (14+ 122/2.8x10"Y?, (15) 
where linearly polarized light has been taken. Typical values of y range from 1.2 to 19. 
The characteristic frequency with which electrons gyrate (see Figure 7) about the 
magnetic field is the so-called electron cyclotron frequency @ .¢, where 
Dg. = By /my (16) 


As shown, @ce becomes a significant fraction of 9 , the laser light frequency. 
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Figure 7. Cyclotron motion of an electron. 


Large self-generated fields can significantly modify the interaction of the laser 
light with the plasma [Ref. 15]. This phenomenon can be explained by considering the 
oscillatory motion of the electrons, which results from the electric field of the incident 
laser light. The Lorentz force associated with this oscillation of electrons across the 
magnetic field is determined by : 

F = -€ os X Bac . (17) 
In this equation, Uo; is the velocity of electron oscillation in the laser field. This Lorentz 
force can become very strong. For Bac= 10° Gauss and taking the peak v,; = c, the force 
corresponds to an oscillating electric field with amplitude of approximately 3x10'' V/cm! 
This large force can be expected to enhance the absorption efficiency, increase the heated 
electron energies, and lead to the generation of harmonics of the incident light. Fora 
discussion of some of these effects at much lower intensities, see Ref. 16. 


To demonstrate the effects of the self-generated fields, we used so called 
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1 1/2 D simulations (1 position, 2 velocities i.e. velocities in two directions for electrons 
and ions) with Zohar, a Particle—In-Cell code (PIC) that solves the complete set 


Maxwell’s equations and uses relativistic dynamics [Ref. 7]. 


As shown in Figure 14, an intense light pulse was sent from vacuum onto a slab of 


overdense plasma. 


Bac y 
Plasma region 
2 
. Bac 

Eos 

Pe vacuum 
cident —~ 

laser light 
wave 360 a)" 


, 


Pulse Shape 
Figure 14. A schematic illustration of a computer simulation of intense laser light 


interacting with an overdense plasma. 


In these simulations, a dc magnetic field pointing in the y direction was added in the 
code. The oscillating electric field of the laser light E,, was chosen to be either parallel 
or perpendicular to Ba, ( by simply rotating the electric vector of the laser light 90° ). 
Some typical parameters are shown in this figure. The initial plasma density rose rapidly 
to 50x ng, where nz, is the critical density. The intensity of incident light rose to its peak 


and then returned to zero in a time of 360@,'. This corresponds to a pulse width of 180 


fs (femtoseconds or 10"'° seconds) for 1.06pm light. 
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To most simply illustrate the effects, consider an example in which the 1.06 um 
laser light intensity is 107° W/cm’, the magnitude of the imposed Bg, is 300 MG, and the 
ions are fixed. The enhanced heating due to the dc magnetic field is apparent in Figure 
(15), which shows electron phase space from simulations with and without the external 
Bac field. More electrons are heated to higher velocity with the B field, as expected due 
to the additional force (-e vs x Bg,). Similar results can be seen from the distribution 


functions (number of electrons vs. energy) shown in Figure (16). 
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Figures 15 a,b. Snapshots of electron phase space in simulations a) with Ba, and b) 


without Bg. 
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Figures 16 a, b. Plots showing numbers of electrons versus energy a) with Ba, field 


on b) without Ba, field. 


Absorption can be measured by several different methods. In our approach, we 
chose to measure absorption by computing the change in kinetic energy of the plasma 
(from times before and after the laser pulse) and then dividing by the energy of the laser 
pulse. Figure 17 shows the absorbed energy versus time. In these simulations, the upper 
curve corresponds to the simulation with the magnetic field and the lower curve is the 
simulation without the magnetic field. If we normalize to the energy of the laser pulse, 
the fractional absorption was approximately 19% without B,,; it rose to about 30% with 


Bac. 
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We also ran a simulation with Bg:, but with it pointing in the Z-direction. In this 
orientation v,; and Ba, are parallel, meaning that vo; x Bac = 0. As we expected, there was 


no enhancement of the absorption observed. 


Bac = 0 


fabs = 19% 
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Figure 17 showing enhanced absorption when Bg, is turned on. 


The simulations also show that the force due to oscillations of electrons across Bac 
leads to an intense generation of even harmonics in the reflected light (see Figure 18). 
These even harmonics arise from the beating of the incident light wave with the density 


fluctuation driven by the v,, x B,. force. In addition, there is enhanced radiation at lower 
frequency i.e. @ <@,, which we attribute to radiation by electrons as they gyrate about 


the strong Ba; field. It would be interesting to look for this radiation in experiments, as 
detection of this radiation may serve as a diagnostic for the large self-generated fields that 


we propose can be generated in the ultra-intense laser regime. 
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Figure 18. Frequency components of the reflected electric field showing even 


harmonics with Ba. 


Lastly, we carry out the same simulations allowing the ions to move. The 
electron-ion mass ratio is taken to be 1/3700 (i.e. m./mj= 1/3700). As shown in Figure 
19, we again find that the absorption is greater with Bg, (= 20%) than without Bac (fabs® 


2%). It is somewhat reduced in both cases because the intense pressure of the incident 
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light wave strongly steepend the density profile of the plasma. Again, we find that the 
heated electrons are significantly more energetic in the simulation with Bac. This is clear 


from the electron phase spaces shown in Figure 20. 
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Figure 20. Plots of electron phase space in simulations with moving ions a) with 


and b) without Bac. 
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The ion phase space shown in Figure 21 illustrates the ion acceleration into 
vacuum on both sides of the plasma. This acceleration is produced by the electric field 
that sets up to prevent the electrons from leaving the ions behind. Note also in Figure 21 
the strong ion acceleration as the pressure of the intense light pushes on the plasma. With 
Bac, ions that are more energetic are generated both in the expanding plasma and at the 


light-plasma interface. 
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Figure 21. Simulations of Ion phase space with and without Bg, 


a2 


Finally, the simulations show a very interesting compression of the Ba. The 
pressure of the intense light wave compresses the density of the plasma, as it reflects (and 
partially absorbs). The peak magnitude of the field can be significantly increased. As 
seen in Figure 22, the peak By, has gone up by a factor of about three, which corresponds 


to over 10° Gauss in this particular example. 
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Figure 22. Simulation of the spatial profile of the magnitude of the self-generated 
magnetic field after irradiation by an intense laser pulse with intensity of 10” 


W/em’. 
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IV. CONCLUSION 


We have investigated via computer simulations how self-generated magnetic 
fields up to 10’ Gauss affect the interaction of short, ultra-intense pulses of laser light 
with plasmas. We have shown that the laser absorption and the high-energy electron 
generation can be significantly enhanced by these self-generated magnetic fields. The 
enhancement is eres to the very large Lorentz forces resulting from the oscillation 
of electrons by the laser light across these fields. 

Other effects shown are the generation of intense harmonics in the reflected light 
and noticeable cyclotron radiation. Finally, the large self-generated magnetic fields can 
be compressed to even large values by the pressure of the ultra-intense laser light. 

Since the late 1970’s, the existence of self-generated magnetic fields, in laser- 
irradiated targets, have been investigated. Some of the initial research confirming the 
existence of these magnetic fields was conducted by scientists at the Naval Research Lab 
[Ref. 17]. To attain magnetic fields on the order of 10° Gauss, which we have 
considered, requires a powerful laser capable of irradiating a target with intensities on 
the order of 10?! W/cm”. Already, such intensities have been achieved by the 
“Petawatt”, short pulse laser at the Lawrence Livermore National Laboratory. 

Exciting high-energy applications and new research in creating Antimatter using 
ultra-intense lasers like Petawatt laser is being conducted. As an example, intense light 
from the Petawatt laser has been directed onto a thin gold film where it creates a plasma 
plume, which acts as a sort of messy Wakefield accelerator [Ref. 18] (see Figure 23). In 
particular, the electric field of the laser light rips electrons from the gold atoms and send 


the electrons shooting off with energies as high as 100 MeV. 
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Figure 23. Illustration of an ultra-intense laser matter interaction. The laser’s 
electric field is strong enough to rip electrons from the gold atoms and shoot them 
off with energies up to 100 MeV, initiating the generation of y rays, which in turn 


can initiate fission or the generation of anti-matter (electron-positron pairs), From 


Ref.[18 ]. 
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The generation of the relatively lower 3-10 MeV electrons using ultra-intense 
laser light (intensity about 107’ W/cm”), is an important issue for laser radiography. 

Some of the higher energy electrons radiate gamma rays, which in turn can induce fission 
or create electron-positron pairs (the first antimatter made in laser-solid interactions). A 
positron is an electron with a positive charge and is the anti-particle of the electron. The 
positron does not exist naturally on our planet. It ‘arrives’ in cosmic radiation or is 
created in the laboratory. Thus, laser photons at the electron-volt level can, by teaming 
up, initiate the sort of million-electron-volt nuclear reactions that normally take place 
only in an accelerator. Moreover, the femtoseconds laser pulses can be focused to a 
much smaller spot size than is possible with any conventional particle beam . 

Continued enhancements and improvements in technology offer rich and exciting 
frontiers of science to be explored. Continued probing and exploring of the ultra-intense 
regime of laser plasma interactions shall give more understanding of very high 
temperature matter. This understanding will define and point to new applications, 
ranging from developments in laser radiography to the successful application of the fast 


ignition concept of inertial confinement fusion. 
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